A NOTE ON THE KINEMATICS 
OF DISLOCATIONS IN CRYSTALS. 



Jeffrey Comer, Ruslan Sharipov 



Abstract. A part of the theory of dislocations in crystals is revised with the aim 
to fit it into the framework of the nonlinear theory of plasticity initially designed for 
amorphous glassy materials. 



1. Geometry of dislocations. 

Dislocations provide a microscopic mechanism explaining the plasticity of crys- 
tals. The idea of dislocations suggested by Taylor and Orowan in 1934 is illustrated 
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Fig. 1.1 



Fig. 1.2 





Fig. 1.3 



Fig. 1.4 



on the first two figures. Fig. 1.1 and Fig. 1.2 show screw-type and edge-type dislo- 
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cations respectively. There are also mixed type dislocations combining the features 
of both screw and edge dislocations. 

In the continuous limit, a single dislocation can be understood as shown on 
Fig. 1.3 and Fig. 1.4 above. Imagine a smooth curve AB within the continuous 
medium being the edge of some surface S. This fact is denoted as AB = dS. 
Imagine that the medium is cut along the surface S and then glued with some 
displacement (see Fig. 1.4). Upon gluing, all points of the surface S outside the 
dislocation line AB become regular points of the medium (as regular as all other 
points within the medium). Though a dislocation produces the stress and clastic 
deformation around itself, it doesn't produce the defects of crystalline grid outside 
the dislocation line AB. Therefore, if we cut out a sufficiently small spherical 
neighborhood G of some point C ^ AB and then release it, we get some stress-free 
crystalline body with no defects of the crystalline grid (see Fig. 1.5 and Fig. 1.6). 
Mathematically this fact is expressed by a map from the spherical neighborhood G 




Fig. 1.5 Fig. 1.6 

of the point C to some domain fl in M 3 (see Fig. 1.6). This map can be given by 
three functions 

x 1 = x 1 (y 1 ,y 2 ,y 3 ), 

x 2 = x 2 (y\y 2 ,y 3 ), (1.1) 
x 3 = x 3 (y 1 ,y 2 ,y 3 ). 

Here y 1 , y 2 , y 3 are some curvilinear coordinates within the real crystalline body, 
while x , x 2 , x 3 are Cartesian coordinates in the three-dimensional space R 3 which 
we used in our mental experiment where we cut out the ball G. The Jacobi matrix 

T ! = % < L2 » 

of the map (1.1) is non-degenerate since otherwise it would mean the infinite com- 
pressibility of the medium. The map / : G — > Q is invertible and due to the 
non-degeneracy condition det T ^ 0, the inverse map f^ 1 : fi — * G is also given by 
three smooth functions similar to (1.1): 

1 1/1 2 T\ 

y = y i {x\x*,x' 3 ), 

y 2 =y 2 (x 1 ,x 2 ,x 3 ), (1.3) 
y 3 = y 3 {x 1 ,x 2 ,x 3 ). 
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By differentiating (1.3) we find the inverse Jacobi matrix S 



en = fly* 

dxi ' 



(1.4) 



Note that the map (1.1) can be extended to any connected and simply con- 
nected domain within the crystalline body that comprises the point C and does not 
comprise the points of dislocation, line AB. Such an extension is non-degenerate 
(detT 7^ 0), but in the general case it is not globally bijective. In such a case, the 
inverse map (1.3) is defined only locally. However, the components of the mutually 
inverse Jacobi matrices (1.2) and (1.4) can be treated as global functions 

Tl=T?{y\y\y% S] = Sj(y\y 2 ,y 3 ) (1.5) 

defined at all points of the crystal except for those lying on the dislocation line AB. 





Fig. 1.7 



Fig. 1.8 



Now let's consider a crystal with one dislocation line AB and define the following 
path integral along some closed path 7 that encircles the dislocation line AB: 



(1.6) 



(see Fig. 1.7). Here r 1 , r 2 r 3 are the components of the tangent vector r of the 
path 7. For the path given parametrically by three functions 

y 1 = y 1 (s), y 2 = y 2 (s), y 3 = y 3 (s) 

this tangent vector is defined by three derivatives 

r> = f, ,- = 1,2,3. 
as 

Note that the path integral (1.6) along any path that does not encircle the 
dislocation line (see Fig. 1.7) is identically zero: 



I fj T J ds = 0. 



(1.7) 
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Indeed, due to (1.2) the integral (1.7) is transformed into the path integral of the 
second kind applied to the total differential of the smooth function x l (y 1 , y 2 , y 3 ): 



Theorem 1.1. The value of the integral (1.6) is an invariant of the dislocation 
AB. It does not depend on a particular contour 7 encircling the dislocation line. 

The proof of this theorem is obvious from Fig. 1.8. Indeed, for the pair of 
contours 71 and 72 on Fig. 1.8 we have 



Definition 1.1. Three constants b 1 , b 2 , b 3 determined by the integral (1.6) are 
the components of a vector b characterizing the dislocation line. This vector is 
called the Burgers vector of a dislocation. 

Note that the Burgers vector b is not a vector in the space of real crystalline 
body. It is associated with the imaginary space of stress-free crystalline matter 
shown on Fig. 1.6. In what follows we shall call this space the Burgers space. 

The concept of Burgers space is convenient for understanding the nature of 
matrices (1.5). Although they have upper and lower indices and depend on the 
coordinates of a point in the real crystalline body, they are not components of 
traditional tensor fields. They are double space tensors. The index j in T? is 
a covariant index associated with the space of real crystalline body, while i is a 
contravariant index associated with the Burgers space. As for the inverse matrix 
S = T _1 in (1.5), its lower index j is associated with the Burgers space, while its 
upper index is related to the real crystalline body. 

Usually each dislocation line is a closed path within crystalline body. Otherwise, 
if it is not closed, it should begin at some point on the boundary of the crystalline 
body and it should end at some other point which is on the boundary as well. 



Usually, each dislocation line is taken with some orientation assigned to it. One 
can change the orientation of a dislocation line, however, in this case its Burgers 






Fig. 1.9 



Fig. 1.10 
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vector b is changed for the opposite one: b — > b. If the closed contour 7 encircles 
several dislocation lines (see Fig. 1.9), then 

I fj r j ds = ±b\l) ± . . . ± b l {N). (1.8) 

i, 3=1 

The sign of each Burgers vector in right hand side of (1.8) is determined by the 
orientation of corresponding dislocation line. 

In some cases dislocation lines have brunching points as shown on Fig. 1.10. For 
the Burgers vectors of dislocation lines in this case we have the equality 



bi=b 2 +b 3 . (1.9) 

The equality (1.9) is an analog of Kirchhoff rule for currents in electromagnetism. 
Its proof is clear from Fig. 1.10. 



2. Continual limit for dislocations. 

In order to detect macroscopic phenomena associated with dislocations we should 
have a substantial amount of dislocations in each macroscopically essential volume 
of the medium. In this case, instead of considering individual dislocation lines, we 
consider the density of Burgers vectors for dislocations 

Pi = pi(y\y 2 ,y 3 ). (2.1) 

Like S and T in (1.5), the functions (2.1) are components of a double space tensorial 
field. The upper index i in (2.1) is associated with the Burgers space, while j is a 
traditional tensorial index associated with the space of real crystalline body. 

Since dislocation lines cannot end in the interior of a crystal and since they obey 
the conservation law (1.9) at their brunching points, the amount of Burgers vectors 
flowing into some domain trough its boundary with dislocation lines is equal to 
the amount of this vector flowing out of this domain. This fact is written as the 
following integral equality for the density of Burgers vectors (2.1): 

3 

Y,pj njdS = °- ( 2 - 2 ) 

on i =1 ' 

Here n , n 3 , n 3 are components of the unit vector n of the external normal to the 
boundary <9£1 of the domain f2. The differential form of the equality (2.2) looks like 

3 3 

EE^' V ^ = °- ( 2 - 3 ) 

J=l fe=l 



The covariant derivative V^p 4 in the formula (2.3) is calculated as follows: 
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Note that g 3 in (2.3) are components of the metric tensor and 1^ ■ in (2.4) are com- 
ponents of the metric connection. They are determined by the choice of curvilinear 



coordinates y 1 , y 2 , y 3 (see [1]). Note also that in writing (2.4), we do not apply the 



standard rule of covariant differentiation to the upper index i. This is because it is 
associated with the Burgers space other than the space of real crystalline body. In 
short form the equality (2.3) is written as 



div p = 0. 



(2.5) 



In this form (2.5) the equality (2.3) resembles the Maxwell equation divH = in 
electromagnetism. However, one should remember the difference: H is a vector, 
while p in (2.5) is a double space tensor. 

Let S be some surface spanned onto the contour 7. Then 7 = OS (see Fig. 2.1). 
The total flow of Burgers vectors across the surface S is given by the following 

surface integral (compare with (2.2)): 




i=i 



(2.6) 



The value of the integral (2.6) cannot 
change unless some dislocation lines move 
and cross the contour 7. For this reason 
the time derivative db % j dt is given by some 
path integral along the contour 7: 



Fig. 2.1 



db l 



(2.7) 



fe=l 



The double space tensorial quantity j£ in (2.7) is called the current of Burgers 
vectors produced by moving dislocations. Combining (2.6) and (2.7), we get 



s j=1 as k=1 



(2.8) 



The integral equation (2.8) can be transformed to differential form by applying the 
Stokes formula. As a result we obtain the following equation: 



q—l p—1 r—1 m— 1 



(2.9) 



Here g qr and g pm are the components of metric tensor, and ujk qp are the compo- 
nents of a completely skew-symmetric tensor. It is called the volume tensor. Its 
components are expressed through Lcvi-Civita symbol: 



^kqp — J dct(<7jj) Ekqp 
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(see more details in [1]). When applying the covariant derivative V r to a double 
space tensor one should remember that the indices associated with the Burgers 
space are ignored. For the derivative Vrj^ in (2.9) we have 

Vrj4 = §^-£r* m j< (2.io) 

" s=l 

(compare (2.10) and (2.4)). Like the equation (2.5) above, the differential equation 
(2.9) can be written in a shorter form: 

|£+rotj = 0. (2.11) 

Like the equation (2.5), this equation (2.11) is an analog of corresponding Maxwell 
equation in electromagnetism (see [2 ]) : 



1 <9H „ „ 
- — + rot E = 0. 

c at 

Now let's return to the matrices (1.5). For a single dislocation they are de- 
termined by formulas (1.2) and (1.4). Usually they arc singular functions at the 
points of dislocation line, just like Coulomb potential of a point charge. However, 
if charges are treated as continuously smeared in the space, the electric potential 
(p(t, y , y 2 , y 3 ) is a smooth function. Similarly, in continuous limit of dislocation 
theory TJj (i, y 1 , y 2 , y 3 ) and Sj(t, y 1 , y 2 , y 3 ) are smooth functions forming two ma- 
trices inverse to each other. In this case they are not determined by formulas (1.2) 
and (1.4) any more. Instead, we have the equality 

<fj2^rUs=[j2p> j dS (2.12) 

OS i = 1 s ^ 

derived from (1.8). Applying the Stokes formula to (2.12), we get 

rotT = p. (2.13) 
In coordinate form the equation (2.13) is written as follows: 

E E E E 9 qr 9 pm V r 2£ = pi 

q—1 p—1 r—1 m— 1 

Comparing (2.13) with H = rot A, we conclude that the double space tensor field 
T with components Tj(t, y 1 , y 2 , y 3 ) here plays the same role as the vector-potential 
A in electromagnetism. 

3. Deformation tensors. 



Let's consider the motion of a crystalline medium in the presence of dislocations 
in it. Here we reproduce in part the content of [3] in order to have the same 
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notations. Suppose that a point of the medium with coordinates y , y 2 , y 3 has 
moved to the point with coordinates y , y 2 , y 3 . Then we have a map: 

y 2 = y 2 (t,y\y 2 ,y 3 ), (3.1) 
y 3 = y 3 (i,y\y 2 ,y 3 ). 

This is the displacement map r. The argument t in (3.1) is responsible for the time 
evolution of the displacement. The time derivatives of the functions (3.1) determine 
the components of the velocity vector v: 

= i = 1,2,3. (3.2) 

As defined in (3.2), v , v 2 v 3 are the functions of t, y , y 2 , y 3 . However, in order 
to interpret them as the components of a vector field, they should depend on the 
coordinates of the current actual position of a point of the medium. To change the 
arguments of the derivatives (3.2) we use the inverse displacement map r _1 : 

y 2 =y 2 (t,y\y 2 ,y 3 ), (3.3) 
y 3 = y 3 (i,2/\y 2 ,y 3 ). 

The time dependent maps (3.1) and (3.4) define two Jacobi matrices S and T: 



dy l 
dp 



fi = d & 
J dyi ' 



(3.4) 



The nonlinear deformation tensor G then is defined as follows (see [3] ) : 

3 3 

G, =Y.Y.9rs{y\y\f)Tlf;. (3.5) 



r— 1 s—1 



Upon transforming all arguments in (3.5) into t, y , y 2 , y 3 we get a tensor field G 
with components Gy = Gij(t, y , y 2 , y 3 ). Differentiating (3.5), by direct calcula- 
tions one can derive the following formula: 



9Gk 

at 



v k V fc Gy ^ G kj V lV k - J2 C - ! V J 



(3.6) 



fc=i 



fc=i 



fe=i 



In order to describe the plasticity of amorphous materials in [3] the following de- 
composition of the deformation tensor G was suggested: 



3 3 

Gij = 2J ^2 Gi G kq G 1 -. 

k=l q=l 



(3.7) 
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Here Gkq are components of the elastic deformation tensor G, while G k and G? 

are components of the plastic deformation tensor G. For these tensor fields G and 
G in [3] the following evolution equations were suggested: 



dG 



kq 



3 



dt 



J2 v r V r G feg = -J2 ^kv r G rq - J2 &kr V,i 

r—1 r—1 r—1 

3 3 

+ ®k Grq + 2J Gkr Qq- 



(3.8) 



dG[ 
dt 



r=l r=l 
3 



Y, v r V r G? = (G[ V r v k - ViV r G r k ) - r G l- ( 3 - 9 ) 



The main goal of the present paper is to show that the decomposition (3.7) and 
the differential equations (3.8) and (3.9) can be consistently incorporated into the 
existing theory of plasticity in crystals in its nonlinear version. 

4. Kinematics of a dislocated medium. 

Let's begin with the equation (2.11) and substitute (2.13) into it. As a result we 
obtain the following differential equation: 




It is known that a vectorial field with zero curl is the gradient of some scalar field: 

— +j = -gradw. (4.1) 

In our case all of the fields T, j and w in (4.1) are double space tensors; they have 
one upper index associated with the Burgers space. Therefore, (4.1) is written as 

-£ + *—w (4 ' 2) 

The vector w in (4.2) can be interpreted as the velocity vector (it easy to check that 
its components w , w 2 , w 3 are measured in cm ■ sec -1 ). However, w 7^ v. Indeed, 
the components of the velocity vector v defined in (3.2) and then transformed into 
the arguments t, y l , y 2 , y 3 by means of the map (3.3) represent a traditional tensor 
field with one upper index, while w is a double space tensor field. Below we shall 
understand w as an independent parameter of a dislocated medium. The physical 
nature of this parameter is not yet clear to us, it will be studied in separate paper. 
However, there is a transparent analogy with electromagnetism: 

- — + E= -gradp. 
c at 
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Comparing (4.1) with this equality, we find that w is an analog of the scalar po- 
tential ip of electromagnetic field. Moreover, this equality supports our previous 
associations of T with A, j with E and p with H. 

Suppose that the initial state of our crystalline medium is free of dislocations. 
Below we assume that it is stress-free too. Then 

pl(0,y\y 2 ,f) = 0, ji(0,y\y 2 ,f)=0. (4.3) 

Due to (4.3) we can arrange the bijective map from the space of the real crystalline 
body to the Burgers space. It is given by three functions 

' x 1 = a; 1 ^ 1 , j/ 2 , j/ 3 ), 

< x 2 = x 2 (y 1 ,y 2 ,y 3 ), (4.4) 
_ x 3 = x 3 (y\y 2 ,y 3 ). 

The map (4.4) is an isometry because we assume that the initial state of the crystal 
has no deformation. The isometry condition is written as 

3 3 

g rs (y\f,f)=Y,Y,9 Pq T p r h (4.5) 

p=l q=\ 

Here *g pq are the components of metric tensor in the Burgers space, while T£ and 
Tf are the components of Jacobi matrix for the map (4.4): 

r dy^ 

Note that g pq = const since we choose Cartesian coordinates in the Burgers space 
(see Fig. 1.6). 

The next step is to add the time variable to the map (4.4). For this purpose let's 
use the inverse evolution map (3.3) and let's consider the composite map 

' x 1 = x 1 {y 1 {t,y 1 ,y 2 ,y 3 ), ... , y 3 {t, y 1 , y 2 , y 3 )), 

< x 2 =x 2 (y 1 {t,y 1 ,y 2 ,y 3 ), ... ,y 3 (t,y\y 2 ,y 3 )), (4.6) 
_ x 3 =x 3 (y 1 {t,y 1 ,y 2 ,y 3 ), ... , y 3 (t, y 1 , y 2 7 y 3 )). 

Using the chain rule, for the Jacobi matrix of the map (4.6) we write 

dr p * ~ 

y r=l 

From (3.5), (4.5), and (4.7) one easily derives 

<■■, Y.Y.ih.irr. (4.8) 

r— 1 s — 1 
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The equality (4.8) means that the deformation tensor G can be defined through the 
composite map (4.6). As for the matrix (4.7), we interpret as the components of 
a double space tensor T. Both tensors T and T are called the distorsion tensors (see 
|[4-7]| ) : T represents the compatible distorsion since it is given by partial derivatives 
in (4.7), while T represents the incompatible distorsion since the equality (1.2) is 
not valid upon passing to the continuous limit (see (2.12) and (2.13), see also the 
comment above the formula (2.12)). 

The compatible distorsion arises due to the macroscopic deformation of a crys- 
tal. In the elastic case, the macroscopic deformation is transferred to the micro- 
scopic level and produces the same distorsion of interatomic bonds (see Fig. 4.1 and 
Fig. 4.2). The plastic deformation is that very case, when some interatomic bonds 

distorsion 




Fig. 4.1 Fig. 4.2 Fig. 4.3 

get torn and then relinked in a different way. On Fig. 4.3 we see the birth of a 
pair of the edge dislocations with mutually opposite Burgers vectors. On Fig. 4.4, 
Fig. 4.5, and Fig. 4.6 one of them moves from the left to the right. Behind the mov- 




Fig. 4.4 Fig. 4.5 Fig. 4.6 

ing dislocation the series of undistorted cells arises, while the total (macroscopic) 
distortion angle remains unchanged. This fact explains why T / T for plastic 
deformations. 

The elastic response of a body is determined by the elongation and / or contraction 
of interatomic bonds within it. For this reason let's define the elastic deformation 
tensor G by analogy with (4.8), but using T instead of T: 

P=l 9 =1 

Thc components of the plastic deformation tensor G arc defined as follows: 



P =i 



(4.10) 



12 JEFFREY COMER, RUSLAN SHARIPOV 

Here are components of the inverse matrix S = T _1 . Though defined through 

the double space tensors, the deformation tensors G and G are traditional tensor 
fields in the space of the real crystalline body. The indices p and q associated with 
the Burgers space in (4.9) and (4.10) both are summation indices. They disappear 
when the sums are evaluated. 

Note that now the decomposition (3.7) follows immediately from the expressions 
(4.8), (4.9), and (4.10). Therefore, it is sufficient to derive the equation (3.9). The 
equation (3.8) then is derived from (3.9) and (3.6) due to the decomposition (3.7). 
As the first step in deriving the differential equation (3.9), we differentiate the 
equality (4.10) with respect to the time variable t: 

p— 1 p— 1 

In order to differentiate the inverse matrix S = T _1 in formula (4.11), we use the 
well-known standard formula S' = -S T' S: 

q—l r—1 

Substituting (4.12) into (4.11) and applying the formula (4.2), now we derive 

^ = E E ^ U + 9 £r) Gl + ± *£. (4.13) 

q=l r =l ^ y ' p=\ 

For the time derivative of Tf in formula (4.13) we have the following equality: 

9t f-[ dy i y ' ' 

The equality (4.14) is derived in few steps by applying the chain rule to the mapping 
functions (4.6) and to their inverse mapping functions 

V 1 = y 1 (t,y 1 (x 1 ,x 2 ,x 3 ), ... , y 3 (x\x 2 ,x 3 )), 
y 2 = y 2 (t,y 1 (x 1 ,x 2 ,x 3 ), ... , y 3 {x l ,x 2 ,x 3 )), 
k V 3 = y 3 (t,y 1 (x 1 ,x 2 ,x 3 ), ... , y 3 (x\x 2 ,x 3 )). 

Apart from (4.14) we need the identity 

OTP dTf 



Qyl QyT 



(4.15) 



which follows immediately from (4.7). Now, applying the formulas (4.14) and (4.15) 
to (4.13), we derive the following equality: 

q=lr=1 T=X (4.16) 

3 3 nrpp 

p=lr=l y 
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The last term in (4.16) can be transformed as follows: 



3 3 n^P 3 3 £)( qk TP "\ 3 3 o ok 

EE^S = EE^^^-EE^?r 

p=l r=l y p=l r=l y p=lr=l y 



(4.17) 



3 Ptr<k _3_ 3 3 nrpq 

= E^fr+EEE^ r S^- 

r=l w g=i p =i r=l tf 
Upon substituting this expression into (4.16) we derive 

dGl r dG k A / - r dv k dv r - k 

Now let's introduce the following notations: 

^ = ^ + E ^ -^+E^J- ( 4 - 18 ) 

Then the differential equation (4.17) for plastic deformation tensor is rewritten as 

d -§ + E « r d -§r = t (5 % - ^ - g ar. (4.i9) 

By direct calculations one can verify that all of the partial derivatives in (4.19) 
can be replaced by covariant derivatives. As a result (4.19) takes the form of the 
equation (3.9), which was derived in [3] for amorphous materials. 

The partial derivatives in formula (4.18) can also be replaced by covariant ones. 
As a result this formula is rewritten as 

6 r k = V r V k Yl ^ Jr - E ^ + E E ^ ^ V P f r- ( 4 ' 2 °) 

q—1 q—1 q—1 p—1 

Indeed, if we remember that 



_ 8v r r V ' r dv r ' P dv p 2^ L r P 1 m 

" p=l " * m=l 

(see (2.4) and (2.10) for comparison), we see that the Christoffel symbols T k p do 
cancel each other when substituting the above expressions into (4.20). 

The formula (4.20) written in terms of covariant derivatives reveals the tensorial 
nature of the quantities 9 k : they are the components of a traditional tensor field 
(not a double space tensor unlike T, p, and j). 
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Theorem 4.1. The nonlinear deformation tensor G in the theory of crystalline 
dislocations admits the decomposition (3.7) into elastic and plastic parts G and G , 
both satisfying the same differential equations (3.8) and (3.9) as in the theory of 
plasticity for amorphous materials. 

5. Conclusions. 

The theorem 4.1 proved by the above calculations is the main result of present 
paper. It is a purely kinematic (i. e. geometric) result. Among other results, one 
should mention the concept of the Burgers space. The interpretation of Burgers 
vectors as the vectors of a separate space (and, hence, the use of double space 
tensors) is our methodical achievement (in linear theory this feature is completely 
hidden, see [8]). The further development of this technique and the further com- 
parison of amorphous and crystalline plasticity theories (including the dynamics 
and thermodynamics of media) will be done in separate papers. 
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